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PREFACE 


The objective of this study was to provide a theoretical basis for evaluating 
correlation trackers to be used with meter class solar telescopes on the 
space shuttle, and to provide a specific evaluation of the expected performance 
of two trackers - a tracker developed by the Bendix Corporation for Marshall 
Space Plight Center, and a representative tracker using CCD arrays. 

The method used was to digitize several high quality granulation pictures 
from Sacramento Peak Observatory, and run computer simulatxons of the tracker 
studied using the granulation data. Performance of the trackers in tracking 
accuracy and in providing telescope focus information was evaluated. In 
addition to specific evaluation of the Marshall and CCD trackers, parametric 
equations were derived from which the performance of any tracker using the 
same algorithms can be estimated. 

Our principal conclusions are that a modified version of the Marshall tracker 
is suitable for tracking solar granulation to the required accuracy, but not 
for providing focus information. Ihe CCD tracker is suitable for both tasks. 
Finally, the correlation tracker must be supplemented by a limb tracker for 
most applications. 

We recommend that development of a CCD tracker specifically designed for solar 
granulation tracking begin, and that development of the Marshall tracker continue . 
The correlation tracker will be an extremely useful addition to solar telescopes 
in space, and has utility in ground-based observations as well. We also recom- 
mend that further study be conducted in the time behavior of solar granulation, 
in order to better define the utility of correlation trackers used without an 
auxiliary limb tracker. 
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1. Introduction 

The purpose of this study is xwo fold. First, to dev elope the general 
theory of correlation tr.ackers used with solar granulation targets, and 
second to investigate the performance of two specific trackers. The two 
trackers investigated are a tracker developed hy the Bendix Corporation for 
ffershall Space Flight Center, and a representative tracker using CCD image 
detectors. The study concerns xrackers as used with meter class telescopes 
in space, specifically with apertures ranging from 65 cm to 125 cm. The 
aspects of tracker operation studied are short-term tracking accuracy when 
the granulation pattern can he considered as fixed, effects of changing 
granulation patterns on tracking accuracy, and methods of monitoring tele- 
scope focus quality by means of the correlation tracker. 

The first phase of the sxudy was spent in digitizing several excellent 
granulation pictures wixh the FDS microdensitometer at Lockheed. These 
pictures were taken at Sacramento Peak using the 75 cm vacuum telescope. 
Pictures taken hy Lockheed and others were studied, alxhough most of the 
data was taken from a picture supplied hy Dr. Dunn of Sacramento Peak. 

The microdensitometer traces were smoothed to remove film grain noise. 
Enhanced data designed to simulaxe space telescope seeing was prepared hy 
converting the xraces to square waves, and then smearing the traces -^Tith the 
point spread function of telescopes of specific apertures. This data base 
was used for calculating the performance of trackers in the fixed granu- 
lation pattern case. 

¥e had hoped to use a granulation movie supplied hy Dr. Dunn to provide 
information on performance of the trackers with time-varying granulation 
patterns. However, seeing effecxs proved impossible to remove, so this 
direct approach was abandoned. Instead, a theoretical approach using 
stratoscope data was used, which we believe provides a valid model of tracker 
performance. 

This study attempted no actual tracker design work other than inves- 
tigation of various scan pattern geometries and analysis algorithms. Ho 
hardware was considered other than as necessary to obtain image detector 
performance parameters. 

Final conclusions include the following; 
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o Q?he Marshall tracker^ preferably vlth some modifications to scan 
geometry and processing algorithms will provide good performance 
in short-term tracking. 

o CCD trackers ■^ri.U provide even better performance in short-term 
tracking . 

o Any Tracker observing only a few granules ’trill shift badly with 
reject to the limb as the granules move. !Ihe present Ifershall 
tracker which observes only one granule is particularly bad in 
this regard. 

o The CCD Tracker provides excellent focus monitoring potential. 

The lyfexshall tracker is marginal. 

o Since ary correlation tracker observing a small part of the sun 

will follow solar rotation,, a means of compensation must be provided 
if stability with respect to the limh is to be achieved. This 
requires a separaTe limb tracker. 

Perhaps the most powerful use of correlation trackers will be to keep 
two separaTe instruments co-aligned, or to operate inside the control loop 
of a separate limb tracker to reduce the image displacemenus from Thruster 
firings or man push -offs. As a sole tracker the correlation tracker is^ 

in general^ insufficient. 

¥e have prepared parametric representations of tracker performance 
■^^hich allow the easy evaluation of the performance level of any tracker 
using the processing methods studied. These were prepared by modeling 
trackers mth various parameters and allowing them to "track" our digiTazed 
granulation data. The Marshall tracker was evaluated in the same way. The 
model' for this tracker is accurate^ except That its 66 image elements are 
assumed to lie in a line instead of a circle. This was done for processing 
convenience, ¥e show in section 2.5 that this has a negligable effect on 
tracker performance. 

Several unanswered questions still exist about tracker performance in 
the time-varying granulation case. ¥e do not know if motions of neighboring 
granules are correlated. Also, the exact appearance of granulation at high 
resolution is unknown, although the high frequency content of our enhanced 
data agrees well with results from speckle interferometry. These limitaTions 
and proposed future studies are discussed at length in the final chapter of 
this report. 
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2. Basic Principles and Design of Correlation Trackers 
2.1. Preliminary Considerations 

Before the performance characteristics of correlation trackers can he 
discussed, we must first define several basic terms and concepts, First, 
a "tracker" is a device which senses the position of a target object within 
the field of vieTr of a telescope^ and automatically adjusts the pointing 
of that telescope to keep the image of the target object stationary in the 
telescope image plane. A tracker cannot be used to bring an object; initially 
outside the field of vdex-r of the telescope within the field of view, therefore 
it must be supplimented by a "target acquisition system" in order to bring 
the target within the field of view of the telescope. This may be accom- 
plished by manually aiming the telescope at the desired target, by using 
an automaxic acquisition system using a second tracker coupled to a wider 
field finder telescope, or by raster scanning the telescope across the sky 
until the target is located, ivbst real acquisition systems combine several 
of these methods - for example scanning the sky in a search pattern until 
a target appears in a finder telescope, locking on to the targex with the 
finder or "coarse" tracker, then finally tracking the target v?ith the main 
telescope once it has been brought within the field of view. 

It should be obvious that, in order to completely auxomate this acqui- 
sition and track sequence, something must he known about the expected 
characteristics of the target, so that the tracker can identify and track 
the target. In general, xhe more that is knoTO about the target, the simpler 
and more affective the trackers can be. If enough is known a-priori about 
the target that it can be differenxiated from any other object at which the 
telescope is likely to point, then in principle target acquisition can he 
completely automatic. In the case of the sun, which we know to be very 
bright, we could simply build the acquisition system so that it locked on 
to any object tdiich produced a photon flux greater than some Threshold 
picked xo eliminate all dimmer objects. Moreover, once the sun is aequired, 
we can take advantage of the fact that the sun is a nearly perfect hrighx 
disc againsx a dark background to construct a limb tracker which operates 
by simply keeping two points on the sharply defined edge of the sun stationary 
in the field of view. 
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In many cases, however, not enough is know about the target’s appearance 
to allow such a dimple system to work. Indeed, all that may be known about 
the target is that its appearance will not vary too rapidly 't/ith time 
after it first appears in the telescope field of view. This is the sit- 
uation in which a correlation tracker is designed to "STOrk. In solar ob- 
servations from space we are forced to use a correlation trackers when the 
field of view of the telescope is not large enough xo allow constant viewing 
of the solar limb, which is, of course, the only solar feature with an 
appearance that can be accurately predicted in advance. The best that can 
be hoped for on a random patch of the solar disc is a sunspot of some 
arbitrary shape, but relatively stable with time. Islost of the time, however, 
we must be content with tracking features found in a granulation pattern 
characteristic of the quiet sun. These patterns change relatively rapidly 
with time (Correlations drop by l/e in about 6 minutes.), but are stable 
enough for most purposes. The only things- we know a-priori about the 
pattern we must track are statistical in nature - we know the EMS inxensity 
fluctuation of the granulaxion pattern, the average brightness of the solar 
disc, and the average spatial dimensions of the individual granules making 
up the pattern. These items are useful for selecting an image detector for 
use with the correlation tracker, but not for any specific information about 
the pattern to be tracked. The exact paxtern to be tracked must be learned 
after the telescope is pointed at the target. 

2.2 Operating Principles of Correlation 
Trackers. 

A correlation tracker works by using the following steps — 
o The tracker is steered onto the target by some independent means . 
o -An image of the properly positioned target is then stored by the 
tracker . 

o Subsequent images of the target are compared to the stored image 
and the telescope pointed so that the live and stored images over- 
lap as precisely as possible. 

o If changes in xhe target appearance with xime make ix difficult to 
detemnne when xhe live and stored images are properly overlapped, -a 
reference image is stored and tracking continues. 
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The errors and limixations , as well as hhe strengths of correlatiion 
trackers are apparent from these operating principles. The principal ad- 
vantage^ of coirrse, is that it is no longer necessary to know exactly what 
the target looks like '(■sden the tracker is designed^ since the tracker has 
the capability of storing a reference image when it is actually pointed at 
the target. The corollary to this is that the tracker cannox store ixs 
reference image until it Is properly pointed ax the target, and thus is 
incapable of target acquisition although it >ri.U continue to track the 
target after initial acquisition has been accomplished by their means. 

It is possible, of course, to pre-load the image memory with a reference 
Image for use in acquisition, but if enough is known about the target to 
permit this a correlation tracker is probably not required. 

Tracking errors arise from several sourses. 

o Random noise in the tracker image sensor or signal processing 
electronics - including imput photon shot noise. 

o Systematic errors - a particular correlation algorithm may show 
two identical Images as precisely overlapped when in fact they 
are misregistered. 

o Errors due to changes in target appearance with time - it is no 
longer possible to find a perfect overlap between the stored and 
live images. 

It follows from the fact that the tracker contains several sourses of eiror 
that the live and stored images are not likely to be precisely registered 
when the update of the stored image occurs. Thus, the new reference image 
will have been stored ax a slightly jiifferent position than the original 
one. As successive updates occur, the aimpoint of the utracker will drift. 
This drift ^d.11 be random walk if the errors are caused by random detecxor 
noise or by isotropic changes in the target with time. The drift can, 
hotrever, tend in a particular direction (for example toward a brighter 
area of the solar disc) if systematic errors are involved. 

2.3. Basic Design of Correlation Trackers, 

Major design considerations for correlation trackers are; 

0 Selection of the image pattern to be sxored to provide xhe best 
compromise between tracking accuracy and memory size. 
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o 


Selecxion of xlie best criterion for •when the t'wo images (stored 
and live) are most precisely overlapped, 
o Selection of an image detecxor to match available light and target 
contrast. 

It should be kept in mind that all of the above are dependent on target 
characteristics^ ■which are usually known only in a sxatistical sense in these 
applications 'where correlation trackers are used. A final design consideration 
is the selection of the optimum interval between updates to minimize aim- 
point shift. In most cases this can be accomplished automatically by up- 
dating whenever the correlation between the stored and live images drops 
below some predetermined ■value. 

Keeping these general design problems in mind, let us examine in detail 
the problem of designing a correlation xracker for use '^Tith a meter class 
telescope in space. First we will list several characteristics of the solar 
granulation targex with ■?hich we will have to work. 

o Low contrast - a few tens of percent difference in intensity is the 
most thax can be expected between the brighxest and darkest features. 
Granules are typically about 2 arc seconds in dimension, 
o A tracking accuracy of .01 to .1 of the typiical dimension of granules 
is- required. 

o * Except for sunspots, there is' generally a lack of long-lived, well 
defined features, 

o Changes in the target pattern Tvith a time scale on the order of 
minutes can be expected. 

o IThere will be no change in aspect angle or image scale with timei 
These conditions are acxually quite easy to accomodate. The chief 
cause for poor performance of correlation trackers is rapid changes in aspect 
angle or image scale, with time (characx eristic of many terrestrial ap- 
plications). These generate large sysxematic errors and the need for rapid 
updating of the reference image, resulting in rapid shift of the aimpoint. 

In the solar case, the only need for updating is due to changes in the 
granizlation pattern i^ith time, and this problem is much less severe, requiring 
an update only once or t'^rice per minute for optimum performance. (The 
selection of the optimum update interval will be described in detail below. ) 
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The high tracking accuracy required composed to typical granule dimensions 
suggests that for tracking times longer than a small fraction of the 
granulation pattern correlation tiine (6 minutes), unacceptable errors will 
possibly occur due to the motions of indi'jidual granules. For accurate 
tracking with respect to the solar limb the tracker image n/ill have to include 
many independent granules so that their, motions will average out to an 
acceptably low value. 

The selection of an image detector ^ri.11 be influenced by the lot;- contrast 
of the granule pattern, the brightness of the imstge, the spatial resolution 
required, and the bandpass of the servo loop with which the tracker is to 
operate. The selection of uhis detector is described in detail in the 
following section. 

2.4 Image Detectors 

The two -cypes of image detectors treated intensively in this study are 
the image dissector and the solid state detector array. 

The solid s-cate detector array makes use of integrated circuit tech- 
nology to place on a single wafer of silicon an array of independent detectors. 
These detector arrays are available in various patterns and sizes. The 
individual detectors are usually read out sequentially hy one of two methods. 

The readout technique used classifies the array as a "CCD" (charge coupled 
device) or a "CID" (charge injected device). 

A third type of solid state detectors in common use is the "silicon 
diode array" or simply "diode array". Tne diode array method of forming 
the detecuors in the silicon wafer is different from the other two, as is 
the readout technique. 

All three detector types have the same basic spectral response and 
sensitivity - that of pure silicon. Good quantum efficiency is obtained 
from 0.4 to 1.1 'jj,m , peaking at about 80^ at 0.8 |j,m. These characteristics 
can, however, be altered by appropriate "doping" of the silicon. It is 
not within the scope of this study to discuss the characteristics of such 
detectors in full generality. Most solid state detectors are capable of 
signal-to-noise ratios (s/n) of about 300:1 to 1000:1 if enough light is 
available. Most detector elements hold between 10-^ and 10' electrons, so 
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their signal-to-noise ratio is limited hy input photon shot noise to about 

1000:1 to 3000:1. Intrinsic readout noise of 1000 electrons is easy to 

obtain 'wixh mosx types ^ and noise as low as 10 electrons has been reported 

for some CGI) devices. Thus, solid state arrays operated near saturation 

are generally capable of performance close to the theoretical limit imposed 

by photon counting staxisxics. A -cypical value for s/N would be about 

1000:1. Typically, individual detector elements can be read out at a 

speed of au least 2 X 10 elements per second, so that an array consisting 

-4 

of 1000 elements could be read out in 5 X 10 ' seconds or less, sufficient 
for a servo bandpass of about 300 HZ. There is no reason iday the array 
cannot be read out at a slower raxe if insufficient photons are available, 
but as will be shown in a moment, this should not be necessary. 

The image dissecxor is a popular detector for use with trackers, largely 

because its scan pattern can be varied at will, giving ix a measure of 

flexibility not available T/ith other image detectors such as solid state 

detectors or television tubes, ihe penalty paid is a substantial decrease 

in sensitivity and signal-xo-noise ratio-- as compared mth other detectors. 

The detector used with Marshall Tracker is an ITT model F4012 with an S-20 

photocathode. The maximum photocathode current density permitted with this 

tube is IQiA/fem _ The approximate signal -to- noise ratio of xhis detector 

at its maximum current densixy varies from 84:1 (guaranteed) to l68:l 

(theoretical saximum) at the elemenx dwell time used by the Lfershall Tracker 

of 1.5 X lO'^ seconds. If the J-Iarshall tracker is used with a different 

-4 

dwell time than 1.5 X 10 seco nds, the above quox ed signal-to-noise ratios 
are multiplied - "by a facxor of 5 X 10 ^ where is the new 

dwell time. For the present tracker %re mil take a typical value of S/U 
to be 100:1. 

Let us now consider the photon flux available for use by the correlation 
trackers^ in order to determine if any limitations are imposed by input 
photon shot noise. The smallest image element that need be considered has 
a linear dimension equal to half the Rayliegh resolution limit for the 
telescope. Tbe use of smaller elements than this would only degrade the 
performance of the xracker, since no new information would he obtained. 

A telescope of aperture D has an angular resolution a of 1.22 X /d 
where X is the wvelength of the light observed. The area covered by a 
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similar resolution element of half this diameter in the image plane is 
2 2 2 

0.292 \ f /d idiere f is the effective focal length of the telescope. How 

2 

if F is the total photon flux in photons/cm /sec striking the telescope 

4 2 I 2 

apemure, the flux in the image plane will he equal to 1.14 X 10 D F/ f 
and the total number of photons/sec falling inxo the half resolution element 

3 2 -5 

is given by 3.34 X 10 \ F . If we take X = 5 X 10 cm and consider the 
photon flux in the vTavelength region from .4 |j,m to .6 [im^ then F = 9.55 X 10”^ 
photons/cm /sec., and the photon flux per half resolution element from this 
wavelength region alone is 7.97 X 10^ photons/sec. Even if we assume that 
only 1^ of these are counted hy the image detector due to telescope optics 
losses and low quantium efficiencies (10^ would he conservative for image 
dissectors and for solid state arrays)^ we still have about 8 X 10^ 

detected photons per second, which would saturate a CCD having a capacity of 

6 -4 / 

10 electrons in 1.25 X! 10 seconds. For a f /D = 50 system the flux at the 

It 2 

image plane is 4.35 X 10 photcns/cm /sec 'roich if only 1^ are counted 

2 

gives an image dissector photocathode current density of JO mk / cm or 7 OOC 
times the ratial ms.-rimum for the F4012 image dissector. An attmuation 
of at least this amount will he required to prevent damage to the detector. 

From these- figures we can see that for apy reasonable sample rate we 
can assume that either the image dissector or the solid state array ’tall he 
operated at its maximum signal-to-noise ratio. Photon flux provides no 
limitations if a reasonable white-light image is available. If a narrow 
hand blocking filter is used in front of the detector, the image disseexor 
should not he affected down to a bandpass of ahoux l/2A. The solid state 
detector, however, may begin xo lose S/H at a bandpass of between 10 and 200A 
depending on the transmission wavelength and efficiency of xhe blocking 
filter. 

2 . 5 . Scan Patxerns 

In general, it is not possible xo perform a correlation test on the 
enxire image accessible to the detector, resolved to the full power of the 
telescope. The reason is that there is usually so much data to he processed 
in this case that the xracker bandwidth is severely restricted hy the pro- 
cessing time required. Also, the memory required xo store such a large 


2-7 



nvunber of data poinxs may be excessively expensive both in dollars and power 
consumption, althou^ both these ^factors are becoming less important as 
large-scale integration technology continues to advance. 

Given xhax all of xhe points in the image cannon be recorded and anal- 
yzed, it becomes necessary to select some subset of image elements which 
gives adeauaxe information for the req.uired tracher accuracy, while at the 
same time keeping the total number of recorded elements small enough to 
allow sufficiently rapid processing of the data. Since only a fevr of the 
total mmiber of available image elements are selected, it becomes important 
to select those portions of the image ■which contain the most information 
about image position. If the exact nature of the image is knoTm before- 
hand this task is easy. Since an image detector measures light intensity 
as a function of position in the focal plane and since the tracker infers a 
change in image position by noting changes in intensity at certain of those 
positions in the focal plane, it is clear that we wish to monitor those 


parts of the image where the intensify is changing most rapidly as a function 

Si 

of position ( ^ is large). Also, there may be certain types of image 
motion, sucu as roxation about a known poin't, 'which can be easily detected 


and measured without extensive calculations if a particular scan pattern is 


used. 


Most tracker scan patterns developed in the pasx have been designed to 
track objects such as aircraft or targets on the ground. A characteristic 
of these kinds of targets is a rapid change in aspect angle and target size 
with time. Since xhese conditions require rapid updating of the stored 
images, pure correlaxion trackers tend to drift excessively. If, however, 
xhe target has sufficiently high contrast with respect to the background, 
a proper selection of scan patterns 'will allow tracking on the target edges 
over a 'wide range of aspect angles and ranges. A popular pattern for this 
kind of application is the "rosexte”, shown in figure 2-1. As used -i-ath a 
solar tracker this pattern could be used to track sunspots or pores. A 
rosette patxern could also be used to track the solar disc itself as a limb 
guider. 

It is not so apparent', however, that ary advantage is gained when the 
rosette is used ’wixh a true correlation tracker. In that case the tracker 
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TYPICAL TRACKER SCAN P AHERNS 


B-PCKSfl EEa. 




ROSETTE 


TRACK TO KEEP AB = EE 
AND CD = GH 

THIS IS A WAY TO LIMB GUIDE WITH A CORRELATION TRACKER - 
WILL WORK WITH A WIDE RANGE OF TARGET DIAMETERS. IT IS 
USEFUL WHEN TRACIONG SHARP-EDGED TARGETS AT VARYING 
RANGES AND MAGNIFICATIONS. 

RASTER SCAN 


STORES ENTIRE TARGET IMAGE AT THE RESOLUTION OP OPTICAL 
SYSTEM. A MAXIMUM POTENTIAL PRECISION AND ACQUISITION 
RANGE IS AVAILABLE, BUT BANDPASS MAY BE LOW DUE TO THE 
TIME REQUIRED TO PROCESS SO MUCH DATA. 

CROSS LINES 

THE TECHNIQUE OF THE CCD TRACKER STUDIED, SACRIFICES 
ACQUISITION RANGE FOR SPEED AND SIMPLICITY. 


OFFSET CIRCLES 

THE MARSHALL TRACKER TECHNIQUE - FAIRLY FAST AND 
SIMPLE - ALSO GIVES SOME IMMEDIATE FOCUS INFORMATION. 


Figure 2-1 


is positioned at random on the solar granulation pattern^ and can no longer 
use the simplified edge tracking routine. The question now hecomes whether 
a rosette is more likely to contain a large number of high ^ elements than 

OA 

are equal number of elements in some other arrangement. 

For correlation trackers used with time varying solar granulation 
patterns as targets^ it is necessary to know the area over which the motions 
of the granules are correlated. A tracker which uses a scan pattern extending 
over only one such area will be subject no large errors due to local motions 
of the granules TvTith respect to the solar limit. Strictly speaking^ this is 
a problem of image scale rather than scan pattern^ but a pattern idiich 
covers a large spatial extent for the number of elements (eg - widely scattered 
points) would tend to be more immune to this problem than a compact pattern 
(eg - a circle) containing the same number and size of elements. 

For use with a target consisting of random solar granulation patterns 
the only important parameters are the nunlter of image elements, the number 
of independent granules sampled, and xo a lessor extent the size of the image 
elements in arc seconds. !The shape of xhe scan pattern is important only 
insofar as it affects the number of individual granules sampled. 

Consider a granulation pattern with intensity l(X, Y), and further 
consider the value of ^ at each point in that pattern. Assume that an image 

OA 

detector samples that pattern at U positions (x. , y, ) i = 1, 2, , . . . , IJ. 

Let D be the expectation value of ( ^ ) oyer the entire image. Let 
< "be the expectation value of i ^ ^ ^ ^ pattern 

P is superimposed ax random on the granulation image. The larger xhe value 

2 ' 2 
of < D >, the bexter the scan pattern. Note, however that since < I) > 
p P 

is the most probable value over many random superpositions and not a value 
determinable from a single superposition, that 

N 


^2 1 

< D > = ^ S 
P ^ i=l 


hi 

hX 


C 


= D 


X 


( 2 - 1 ) 


Thus, regardless of the pattern of elements chosen, the probable effectiveness 
of xhat pattern in providing maximum sensitivity in converting image motions 
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to intensity variations is the same^ regardless of the nuniber of elements 
sampled or -cheir geometric relationship. Gdiis result applies to any scene 
■where the oracker scan pattern is placed at random on the scene '^Tithout any 
attempt to locate it in a particularly advantageous spot. 

The above result does not mean that the number of elements and scan 
patterns have no relationship to -cracker performance. A tracker -vd-th a 
large nuniber of elemen-cs ’trill always be more precise than one with a smaller 
number. This is true simply because each element acts as an independent 
tracker. Since the light level at each element is measured only to the 
precision permitted by system noise^ and since the position error signal 
from an N element tracker is the average of all N elements, the probable 
improvement in tracker signal-to -noise is a factor of /"N. 

The scan pattern also has a definite effect on tracker performance. 

The most probable accuracy in tracking on a parcicular stored piece of 


granulation image is independent of scan pattern, but the degree by which 
individual stored images give accuracies departing from the most probable 
depends on the "compactness" of the pantems. That is, if 


^ 9X ‘xi^ 


xhen 





< A> 


2 


(2-3) 


is pattern dependent . 

This result be seen in-cuirively by considering^ the exnreme case where 

all U elemenxs overlap, consisting of Jf independent evaluations of I at a 

particular point. Then if xhat poinx happened to fall on a part of the image 

■where ^ o, the tracker would be very poor. If, on the other hand 

oA oA 

at that point happened to be large, the performance TOuld be excellent. 

A ■wide valuation in performance would be expected from stored image to stored 
image* On the other hand, if the IT elements were 'widely separated compared 


to the image structure scale, the probability that the average -value of 
2 

( ^ ) for all the points would vary much from that for the image as a ■whole 

would be small. Reasonably consisxent performance would be expected. The 

2 

result is that a large is expected for patterns small compared to the 
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image structure^ and a small for patterns large compared to image structure. 
This effect has "been investigated empirically for solar granulation as a 
part of this study. The magnitude of the effect depends on the characteristics 
of the granulation pattern and the size of the imeige elements in arc seconds. 

¥e must also consider the effects of different scan pat xerns as they 
relate to acquisition range - the distance hy which the tracker can he dis- 
placed off target without losing its lock on the target. The largest ac- 
quisition range is available^ of course, with a complete raster scan. In 
this case, the lock can, in principle, he re -established after a displacement 
anywhere within the field of view of the tracker. Crossed lines have a 
relatively long acquisition range in the directions along the two lines, 
provided that the motion does not scan along both lines simultaneously. 
Displacement of more than the order of a granule diameter from boxh lines 
simultaneously will result in a loss of lock. This is because no information 
is recorded concerning the character of the image off the two scan lines. 

The offset circle pattern has the smallest potential acquisition range of 
the trackers considered in this sxudy - equal to about one granule diameter 
in any direction. 

None of the proposed trackers have xhe ability to sense image rotation 
as such. This is not required since the rotational pointing accuracy of the 
IK is currently specified as ± 10 irto 3a accuracy. This gives a maximum 
linear displacement of .1 rr if the roll is about one limb of the sun while 
the correlation tracker is monitoring the opposite limb 1920 tt a'vjay. In 
this case the motion in the field of view of the tracker -vrould appear as 
primarily displacement and could be corrected as such. The differential 
displacement due to rotation across the 6o" field of view of the tracker would 
be a maximum of .006 rf and thus is negligable. 

Finally, a few remarks concerning tracker versatility are in order. 

The above consideration concerning scan patterns apply only when the pattern 
is selected before hand and placed on the im.age at random. The image 
dissector used in the Marshall tracker has the advantage that the size, 
shape and position of xhe scan pattern can be easily varied. An operator 
monitoring a television with a representation of the pattern superimposed 
on the television picture could manually alter the pattern size, position. 
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and shape to track a particular feature. This capahility is not present in 
the current Marshall tracker, hut could he easily implemented. 

In summary, the relevant parameters in scan pattern selecrion are 
element size, number, and degree of compactness of the pattern. The 
geometric shape of rhe pattern is relevant only as in affects the compactness. 
¥e have modeled all trackers with linear arrangements of elements and con- 
sidered only the compactness effects of the circular scan pattern of the 
Marshall tracker. The detailed effects of this approximation, as well as 
effects of changing image scale, number of elements, and element size is 
discussed in detail in chapter 3. 


2,6. Analy s i s Algor i t hms 

Once an appropriate scan pattern has heen selected, and a reference 
image recorded, a method must he developed ’tdiich compares the stored and 
live images and derives from them the required information - how far and in 
what direction the telescope pointing has drifted since the reference image 
was fetored. The selection of the proper algorithm to provide the pointing 
information given the raw image data is potentially the mosu complex design 
problem of all. In practice, however, a few algorithms are commonly used, 
and give satisfactory results. 

In addition to the basic algorithm for deriving uhe pointing information 
from rhe image data, it is possible to pre-process the ineige before providing 
it to the pointing algorithm - for example by Fourier transforming it, filtering 
to provide only high spatial frequency data, and re -transforming to image 
space. Potentially the number of pre-processing schemes is infinite, bux 
again, only a few are in common use. Except for smoothing to rejBove grain 
noise and enhancement routines to simulate space telescope performance, we 
have done no pre-processing in this study. In our model of the Marshall 
tracker, however, Me did separaxe our data into 8 discrete levels in the 
same way as the actual tracker, which does constitute pre-processing of xhe 
image. 

The pointing algorithm used by the Ifershall tracker is one of xhe mosx 
common. It consists of forming the correlation function 




C (6) = S 

i=l 


(X, + 6 , 5Tj) 


L (x. , 
1 ' 



(2-4) 
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where the suiumation is over the IT elements of the detector^ represents the 
stored image, and the live image. In the absence of image motion, C(5) 
iri.ll be a maximum ‘tfhen 6 = 0, provided and are substantially the same - 
that is - provided the granulation pattern has not changed very much since 
ihe time was stored. If the xelescope aimpoint changes between the time 
I was stored and the time the correlation is made with the live image I., , 
then the maximum value of C(6) willoccur at some value of 6 ^ 0. The sign 
and magnitude of 6 for which C(6) is largest determines the direction and 
magnitude of the image shift. 

Unfortunately, this simple correlation function is subject to a number 
of systematic errors. The maximum of C(6') iJiay occur at a position where 
I^ and I^ are not exactly overlapped, even if I^ and I^ are identical. To 
see how this can occur, consider the 1-dimensional integral representation 
of the correlation function: 

2a 


(8) = J I (x) I (x + a) dx 


-2a 


( 2 - 5 ) 


Figure 2-2 shows the results of using this function on various distributions 
I.. = I = I. In case I the maximum of C (a) occurs at S = 0. In case II, 
a net slope in the intensity function extending beyond the sampled area 
gives a maximum at 6 = a. A tracker using this algorithm would believe that 
I and I^ were in registration at 6 = a, when in fact they are not. Case 
III shows a typical granula,tion trace and the resultant correlation peak. 

If there is a net slope from left xo ri^t the peak will be slightly offset 
from 8 = 0 as shown. These systematic errors depend on the nature of the 
granulation pattern, and are discussed at length in chapter 3 . 

An algorithm ^jfaich exhibits no systematic errors of this type is .the 
mean square residual function: 

IT P 

c (a) = s [I (x. + a> y.) - iq yq)] (2-6) 

Si X -LiX 

X=1 

we have evaluated the performance of the Marshall tracker and the CCI) xracker 
using this algorithm. 

Given that some method is used to generate as correlation peak of xhe 
type shown in figure 2-2, some means of identifying the peak is necessary. 
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ST“S 


SYSTEMATIC ERRORS IN THE BASIC CORRELATION FUNCTION 


LOCECMEED 




2a 

C(6) 0! / I(x) I(x + 5) dx 
*^-2a 

IN CASE I, A SINGLE JTEATURE ENTIRE LY 
WITHIN THE SAMPLE AREA GIVES AN 
AUTOCORRELATION PEAK! AT a = 0 
AS IT SHOULD. 



C(6) 



IN CASE n, A FUNCTION WITH A SLOPE 
EXTENDING OUTSIDE THE SAMPLE AREA 
GIVES AN AUTOCORRLEATTON PEAK AT 
a = a, NOT THE REGISTERED POSITION. 




IN CASE m, A TYPICAL SOLAR GRANULA- 
TION TRACE GIVES AN AUTOCORRELATION 
PEAK AT A SLIGHTLY OFFSET POSITION. 
THE lARGER THE AREA SAMPLED, THE 
SMALLER THE ERROR. 


Figure 2-2 


The simplest is to simply find the value of § where C( §) is largest. The 
disadvantage of this approach is that in a flat peak, system noise may 
cause large position measurement errors. Another method is to locate two 
points on either side of the peak where 0(5) is the same, and assume the 
peak to he half-TTay in between. This uechniq.ue is less noise sensitive, 
since the slope of the peak is greater on the sides, and an error in the 
value of C(6) causes a smaller error in 6 itself, hut the method is subject 
to errors caused by asymmetrical peaks. These technigues are illustrated 
in Figure 2-3. The Pfeirshall tracker uses the 2nd method. We have analyzed 
the performance of both trackers using both techniques. 

2.7. Effects of the Tracker Servo Loop 

In this sxudy we have concentrated on the expected accuracy of the 
correlation tracker itself, without reference to the searvo loop with which 
it is operated. It is necessary, however, to briefly consider the per- 
formance of the xrackers as componenxs of servo systems. 

The sources of error in a correlation tracker fall into three classes - 
random noise, systematic errors as discussed in section 2.6, and errors 
due to changes in the granulation pattern with time. The first of these, 
random noise, is influenced by the bandwidth of xhe servo loop ^Tith which 
the tracker is operated. 

Consider the Marshall tracker as an example. It makes one correlation 
measurement in about .01 seconds, with a certain random noise component due 
to amplifier noise, input photon shot noise, and other random noise sources. 

If it is used with a servo with a bandpass of lOHZ, it has the opportunity 
to make several independent measurements in the servo response time. Thus, 
the error contribution from this noise component 'VTill be somewhat reduced 
by an averaging process. 

The general rule is that the random noise error will be reduced by an 
/2f 

amount V— -S where f is the servo 3db bandwidth and f is the sampling rate 
fr s r 

in Hertz. 

The maximum bandwidth for which sampling servo system will easily operate 
is tj'^i so under these conditions we expect the random noise contribution 

^ / 2 - 

to be reduced by a factor of v— — or 0 .^ 6 . The present configuration of 
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LOCATION 


LOCKHEED 


PEAIC VALUE - SIMPLY LOCATE THE 
VALUE OF 6 WHERE C(6) IS MAXIMUM. 
NOISE-SENSITIVE ON A FLAT PEAK. 


5 


TWO-POINT BALANCE - LOCATE TWO 
POINTS ON OPPOSITE SIDES OF THE PEAK 
SUCH THAT C(6) IS THE SAME FOR BOTH. 
■ THE PEAK IS ASSUMED TO BE HALF WAY 
BETWEEN. SUSCEPTIBLE TO PEAK 
ASYMMETRY. 

6 


Figure 2-3 




the Marshall tracker has a servo bandpass of 10 HZ_, so a reduction in random 
noise hy a factor of 0 .i {-5 would "he e 3 qpected. 

Unfortunately, this noise component is not usually dominant, hut servo 
bandpass must he considered for precise calculations. The relationship of 
this noise source to others is discussed exhaustively in chapter 3. The 
results given for the random noise and total tracking errors contain no 
allowance for tracker bandwidth. This allowance must he made as provided 
above '^daen the bandwidth is known. If an optimum servo system is assumed 
the random noise figures should be reduced by a factor of O .56 from those 
given in this study. 

2.8. Focus Information 

The two techniques of deriving focus information from a correlation 
tracker are shown in figure 2-4-. The first is simply to measure the sharp- 
ness of the correlation peak by taking the values of C (8 ) at the center of 
the peak and at uwo points displaced an equal and fixed distance on eiuher 
side. The difference between the peak value and the average of the two 
side values gives a measure of the peak sharpness and thus the -quality of the 
focus. The greater nhe difference the sharper the peak and the better the 
focus, ihis meuhod is easily Implemented from data available from the 
Marshall tracker. Unfortunately it is too inaccurate to use. Also, the 
best possible sharpness of the peak, is scene dependenb , thus this method 
gives only a realtive measure of focus, not an absolute one. 

The second method is to take the Fourier power transform of a line of 
video from the tracker and measure the power at high spatial frequencies. 

The more hi^ frequency power the better the focus. This method is also 
scene dependent, but it is much more sensitive. Used with the Marshall tracker 
ix is marginal, but with the CCD tracker it gives excellent res^llxs. It 
requires, of course that the tracker scan a continuous line of elements to 
provide the proper signal for processing. Isolated elements cannot be used 
in this application. 

Results and theory of focus information derivation are discussed in 
detail in chapter 4. 
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X 


INFORMATION 


LOCKMEED 


CORRELATION PEAK SHARPNESS - 
AVAILABLE FROM PRESENT MARSHALL TRACKER 



Figure 2 ^^)- 


2.9. ®ie M^shall Tracker 

Figure 2-5 shows the operating principles of the Ifershall tracker. 

This tracker uses an ITT /4-012 image dissector with a .001” diameter aperture 
and an S-20 photocathode. The scan pattern is a circle containing 66 elements 
which are digitized into either 8 or 2 levels for analysis by a digital 
processor. It' s simple rate is 98.485 HZ and its servo bandwidth is 10 HZ, 
althoiigh there is no fundamental reason why the sample rate and servo band- 
width cannot be increased by at least an order of magnitude if required. 

It's present demonstrated tracking rate is 7 per second, and this could 
easily "be increased to 100 per second with very little trouble. There 
are plenty of photons no do this (see section 2.4). 

The major drawback of the image dissector is a limited signal-to -noise 
raxio. The major strength is its flexibility in providing easily changed 
scan patterns. The present design does not take advantage of this flex- 
ibility hut as we have seen in section 2.5, so long as it is used as a 
randomly positioned tracker this is xinimportant. This type df detector 
would, however, lend ixself xo a manually adjustable sunspot tracker. 

The most important result of the limited signal-to -noise ratio of the 
image dissector is relatively poor perfoarmance as a focus monitor. The most 
serious drawback of the present tracker design is the size.of the areas 
scanned - a circle 1 " in diameter. This makes the tracker extremely 
susceptible to motions of individual granule. , By simply changing the image 
scale so that the tracker covers a larger area this can be remedied, and 
tracking performance much improved. Performance can also be improved by 
using the least squares residuals algorithm instead of a simple correlation 
funcxion. 

2.10, The CCD Tracker 

The t3rpical CCD tracker we chose to sxudy was a tracker using a pair 
of crossed linear detector element arrays each containing ^ 12 . diodes. Arrays 
such as this are available from every major manufacturer of these devices. 

¥e assumed that each detector element had a s/N ratio of 1000:1, and that it 
was digitised with 10 bit accuracy (1024 levels). The tracker was assumed 
to operate by storixig a single scan of the two lines in a 1024 X 10 bit memory 
to serve as the reference image. The equivalent of the Marshal tracker 


2-20 



2-21 



1, STORES 66 ELEMENTS 
AROUND A CIRCLE 

I (i) i = 1,66 

c * 


3. STEPS DOWN, LEFT 

AND RIGHT, AND 

STORES 66 ELEMENTS 

FOR EACH POSITION 

OF THE CmCLE - 

5 POSITIONS IN ALL: 

I , I , L , Ift , I 
c* u* d* J? ’ r 

66 

i=l 

66 

i=l 

Ay(t) = K(C^(t) - C^(t)) 
SIMILARY FOR Ax(t) 



2. STEPS IMAGE DISSECTOR 

SCAN AND STORES ANOTHER 
66 ELEMENTS 
yi) i = 1,66 


4. MOVES IMAGE DISSECTOR 
SCAN BACK TO CENTER 
AND CONTINUES TO READ 
INTENSITY AT ALL 66 ELEMENTS 
AROUND THE CIRCLE 


6. QUANTIZES INTENSITIES 

INTO EITHER 8 OS 2 LEVELS. 

UPDATES WHEN C (t) £ 0.8 
c 


Figure 2-5 


C (x) •would "be fomed "by 
u '' 

511 

C^(t)=S I (i) (-t-, i + 1) (2-7) 

i=l 

N. 

•where I is the stored image and xhe live image at time -ft-, i denotes 

the i’th diode in the vertical line counting from the hottom. Fe have 
examined this tracker using both the least squares residuals and simple 
correlation func-faion algorithms, and have examined both the peak value and 
t-wo point balance methods of peak location. ¥e have also tested the per- 
fomance of the tracker •when the number and size of the image elements are 
varied. 

The tracker gives excellent results both as a tracker and for focus 
information. The chief drawback is the CCD’s rigid scan geometry. It is 
not well suited for manual pattern control. In general, the performance 
level is much better xhan the Marshall tracker, largely due to the greater 
number of elements and better signal-to-noise ratio. 

Fe have made no attempt in this study xo produce a physical design of 
this tracker, although this could easily be done •with some additional efforx. 
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3. Detailed Tracker Simulation 
3.1. Introduction 

An extensive set of computer simulations of tracker performance has 
"been done to evaluate the noise, systematic, and time-dependent errors of 
three different tracker concepts. In so far as possible, the goal has been 
to perform the same operations on digitized granulation scenes as the 
tracker would in real time. The major problem is the lack of relevant 
obsezvations on the temporal behaviour of small-scale features in the 
solar granulation (as distinct from correlations or large-scale average 
behaviour). Even so, tracker performance on fixed scenes has been studied 
in great detail, and a physically plausible model for the growth of tracker 
errors in time has been developed. In this way, upper and lower •limits to 
tracker errors have been evaluated as functions of tracker design, and the 
ranges implied by these limits indicate the uncertainties present. 

For computer simulation, a tracker is defined by the following properties: 
o Tracking Algorithm: the method of inferring image displacement from the 

real-time measured intensities; three were studied, corresponding to the 
Marshall Tracker (MT), a modified Marshall Tracker (MMT), and the least- 
squares algorithm of a CCD Tracker; 

o SN: signal-to -noise ratio of the detector and associated electronics, 

before 'any quantization; 
o number of detector elements; 

o d: detector element size in arcseconds, which is used as a scale factor 

in image space; that is, when d is varied, the geometric arrangement 
of detector elements in the focal plane is constant, but the effective 
focal length and thus the image size changes; 
o T^: average update time, the time interval between scene storages in 

memory. 

Mote that the geometric pattern of detector elements is not varied: 
linear arrays were always used for their great computational advantages. 

As discussed above, tracker geometry has little effect on the errors calcu- 
lated here. ®ie transfer function of the tracker servo-mechanism has also 
not been specified, although it mil affect the noise errors (see Section 2.7). 
This was done because: (a) intelligent servo design requires detailed know- 
ledge of the spectrum of image displacements which the tracker is to suppress; 
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(■b) comparison of different trackers on their own merits as error sensors 
was possible. 

The first step in analyzing a Tracker was to measure its noise and 
systematic errors when wracking on a fixed (time-independenx) scene. 

Since the solar granulation is a random process, errors were computed for 
each of 100 independent scenes (Section 3.2) and ensemble averages were 
found. Several data sets with varying degrees of image enhancement were 
used to study the effects of telescope resolution and seeing conditions. 

For each scene, tracker detector outputs were computed as a function of 
image displacement; then the tracker algoriwhm allowed calculation of noise 
error (standard deviation) and systematic error (constant for a given 
scene). Repetition of these steps provided the two figures of merit for 
fixed scene tracking: mean noise error and root -mean- square systematic 

error Cg. If whe correlation tracker is intended to keep two different 
instruments (spectrographs, for example) aimed at the same point in identical 
scenes, then and Og are the only relevant errors. 

More generally, though, the accumulation of tracking error in time 
is desired. The model presented in section 3.3 pennits evaluation of errors 
as functions of time or update time, given the _algorithm, Ojj, Og and some 
additional information on the tracker field-of-view. If the tracker is used 
to follow a loEg-llved feature (a dark pore or bright granule, for example), 
then only noise and systematic errors must be considered. However, for 
tracking in absolute coordinates with respect to the solar limb, ithen the 
more difficult problem of drift due to transverse motion of granules in the 
tracker FOV must be attacked. A reasonable model, providing upper and 
lower limits to this drift rate, is presented. 

Specific results for three prototype trackers are presented in Sections 
3 .^ - 3 . 6 . In each case, the errors of any tracker using the same algorithm 
can easily be calculated from the formulae and graphs provided. Section 3.7 
gives a fully -worked example of such calculations, and the last section 
contains a brief summary of correlation tracker performance. 

3.2 lata for Tracker Simulation 

The one-dimensional scenes used in deriving tracker errors were digitized 
from a -sdiite light granulation picture, taken in excellent seeing at the 30 
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inch (jS cm) vacuum telescope at Sacramento Peak Observatory by Dr. E. B, 

Dunn. Several other granulation pictures were digitized and gave results 
similar to the SPO frame^ but they were rejected for extensive calculations 
because of poorer resolution. One hundred lines were digitized, separated 
by intervals of 1.5 ", using a O.O 9 I " square slit on the Lockheed PDS 
lOlOG Microdensitometer. Efech line consisted of 10,000 points at O.OO 76 ” 
intervals. After mild smoothing to reduce grain noise, optical densities 
were converted to arbitrary intensity unins using y = ^■ (results are 
insensitive to the precise value). 

Figure 3-1 shows a sample scene. Ihe auto- correlation functions derived 
from this scene by the model Marshall Tracker (MT) and CCD Tracker are shown 
in Figure 3-2. The CCD Tracker' s FOV lies between the dotted lines on 
Figure 3-1, whereas the MT is one-third as long; these trackers are discussed 
in more detail in Sections 3.i|- and 3,6. Mean-squared residual functions for 
the same trackers and data are shoTra in Figure 3-3. 

The excellent spatial resolution of these observations is shown by the 
fact that the correlation peaks in Figure 3-2 are 2C^ narrower than those 
published in Reference 1, obtained with the Stratoscope balloon-born telescope. 
Since tracker performance virtually always improves when fine details (high 
spatial frequencies) are added to the scene, it is imporrant to estimate 
how much detail has been lost because of finite telescope resolution and 
atmospheric seeing. Although the SPO data is nearly diffraction limited 
(Rayleigh limit = 0 .I 7 "), some degradation from seeing is present, and so 
Tracker errors derived from ground-based data are worst cases for the per- 
formance with a 65 cm or 125 cm aperture space telescope. 

To provide lower limits to nhe expecced errors, the ground-based data 
was "enhanced" by quantizing it to a two-level Square wave, representing 
bright granules and dark lanes with perfectly sharp edges. This was merely ' 
a convenient vray of generating random scenes with approximately the correct 
amount of high frequency detail. Figure 3-4 shows the spatial power spectrum 
of this enhanced daxa along with the recent measurements by speckle inter- 
ferometry of the same one-dimensional power spectrum by Kinahan (Ref. 2; 
consistenx resulus also in Ref. 3). After quantizing, the scenes were folded 
with the point-spread function of a 65 cm or 125 cm telescope before use in - 
tracker simulations; the following one-dimensional PSF was used: 
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and D is the diameter^ 65 cm or 125 c™* ligure 3-4 shot^s that the 65 cm 
and 125 cm -enhanced data' yield plausible lo’wer to '-orac'ker errors. 

Finally^ note that the intensity uni'cs used are arbitrary. All trackers 
were assumed to subtract the average level before processing^ and so the 
only important q^uantity is; 


A I 


rms 


< (I - 
I 


( 3 -- 2 ) 


av 

The effective signal-to-noise ratio for detecting intensity fluctuations is 

SN »- A I 1 and all noise errors depend on this product. For the solar 
rms 

granulation, recent measurements seem to be converging on A 10^ 

disk center (Refs. 1-4) j it increases to 15-2C^ at a heliocentric angle of 
5CPand then decreases, falling to zero within 5 " of the limb (Ref. 5). 


3.3 Time-Dependence of Tracker Errors 

The granulation pattern within the tracker FOV changes in time so that 
correlaxion '5'iath a stored scene decreases. This causes noise errors to 
increase with time and also necessitates updating, storing the new scene 
in memory, t'jhich permanently stores the tracker error at the instan'u of updating. 
Furthermore, motion of features in the tracker POF can cause the aimpoint 
or stable point of the tracker to move -with a varying drift velocity (t). 

The auto -correlation function of the granular intensity pattern has 
been measured by Bahng and Sch’t/arzs child (Ref. l); 

< I (t) I (o)r-> = < I (of> (3-3) 

To = 6.27 minutes; 

where the brackets represent ensemble averages. This has been confirmed by 
recent unpublished work of J. P. Mehltretter. The spatial correlation 
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function^ < I (x, t) I (ojo) >;, broadens slightly as t increases. Biis 
broadening is consistent ^rith transverse granule velocities Vq^ of 1 
km/sec or 0.08 arcsec/min. , root-mean-square, A tracker with large FOV 
containing many independent granules can be analyzed very ac(;urately using 
these observed correlation functions. However^ the present Marshall Tracker 
(Ref. 6 ) sees ax most two distinct granules, and so sizable fluctuations 
about this average behaviour must be expected. It is not possible to 
describe these fluctuations without more specific knowledge about the time- 
dependence of individual small granules, and so this report is limited to 
study of average update times tu and errors. 

If the correlation tracker is used to follow a specific feature on the 
sun, then tracker drift between updates is presumably caused by motion of 
the feature. In this case, only the accumulation of noise and systematic 
errors need be considered. The noise error, (At), must be known as a 
function of the time intearval since the last update, At; explicit expressions 
are derived for each tracker algorithm studied below. Whenever an image is 
stored in memory, a systematic error (rms value 03 ) is added, and, for each 
update after initial acquisition of the feature, a noise error (rms value 

)) is also stored. Therefor^, after Tracking for a xotal time interval 
t = nr^ + At, the ” f eature-follomng " error (t) is given-by 


< (t) > = n [ (t^) ] 


+ 0 ^+ 0 ^ (At), 
s n 


(3-4) 


An optimum update time can al>ra,ys be found to minimize < (t) >, several 

examples of which are discussed below. In actual operation, the tracker should 

update when the correlation peak has decreased by a factor EXP ( - t /t ), 

u. o 

being xhe optimum value, and a distribution of actual update times will 
result. 

If tracking in absolute coordinates with respect to the solar limb is 
desired, then tracker drift must be added to < >. To estimate the drift 


rate, suppose that m Independently moving features are presenx in the tracker 
FOV and xhat the drift velocity (t) is the resultant of xhe m random 
velocities. Then (t) has a mean-squared value of and, in absence of 
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any observational data, the auto-correlation function can he assumed to he 


< (it) (o) > - (3-5) 

¥ith this, it is easy to show that; -che men-squared drift after a time 
interval At is: 


< 1 


(3-6) 



[ At-r^ (1 - e“^'^/'"o)]. 


This model can provide reasonable bounds to tracker drift as follows. 
For an upper limi-c, assume that updating has no effect on the drifi; velocity. 
After a time t = nt^ + At, the total tracker error in absolute coordinates 
E(t) obeys 


< (t) > < < E^ (t) > 


+ 


2 

2Y T 
o o 


[t - (1- 


-t/l 


(3-7) 


>) ]. 


For a lower limit, assume that updating resets V (t) to a new, random value, 
^ich leads to the expression 


< (t) > 




(3-8) 


2nV ^ T 
o o 

m 


u 


(1 - e"'^u/'^o) ] 


T / 

+ __2 o [ At - T (1 - e ^^/^o) ]. 

The previous two equations can be used to bound tracker errors as a function 
of time and to find optimum update times. Clearly, the choice of m, the 
nvimber of independently moving features seen by the tracker, is critical 
in determining the drift rate. In the following analyses, one feature for 
every 1.5 - 2 " in tracker length was assumed, corresponding to a typical 
graniile diameter. However, there is some indication (J. P. Mehltretter, 
private communication) that clusters of several granules may move together. 
If true, this would necessitate a larger tracker FOV for a given acceptable 
drift rate. Two-dimensional simulation using a seeing-free granulation 
movie is necessary to settle the matter. Finally, the effects of solar 
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rotation would need to “be compensated if a constant aimpoint witli respect 
to the solar limb is required for more than about 30 seconds. (Assuming 
a shift of more than 0.1 " is unacceptable). 


3.^1- Marshall Tracker (MT) Simulation 
(a) Mathematical Description 

The tracking algorithm of the MT has been fully described above: see 

Section 2-9* It computes the "up" and "doW correlations, and 
according to 

^ 1 t 

C (y«) = I, (y^i d) I. (yO. (3-8) 

d 

Xll 

Here, I.(y ± d) are the outputs of the i detector element in the stored 

3. O f 

up/down pictures, which correspond to the initial aimpoint y^; (y* ) 

is the real-bime result, corresponding to the real-time aimpoint y^ . 

Tracking error is zero when y’ = y^^* These intensities are quantized in ei^T; 
levels after subtraction of the average. 

The stable point of the MT, y^ , is the value of y' where, in absence 
of noise, 

(y* = y^) = (y’ = y^). (3-9) 

As* 3 numetries in the observed correlation peaks cause the systematic error, 
given by 

where the brackets represent averaging over an ensemble of scenes. 

Waen the tracker is aimed at its stable point;, noise in the measured 
intensities causes it to produce a noisy error signal, Ay^. This error 


signal is 


^u - ^d 

3 C (y' ) / 3 y' 


(3-11) 


Its rms value, a^, is given by 

^2 _ ^ ^ (3-12) 

' 1 B C / 3 y‘ 1 ^ " 
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2 

where cr (c) is the variance of C or due to noise. 

u a 

Appendix A that 


(C) = p [2 2 I ^ Jj , p AI^ ]j 
i X el 


P 


0.698 I 

av 

SN . Al 


In is shown in 
(3-13) 
(3-14) 


Al is the quantization unit^ approximately equal to 5 1^^^ AI^^/8, and p 

is the prohahility of a noise-induced quantization error, 

The preceding equations accurately describe tracker errors in one 

dimension. Since the x and y directions are handled independently by the 

MT, the simulation values of o and cr reported in the next section have 

ns 

been multiplied by . 

(b) Fixed-Scene Results 

Noise and systematic errors were evaluated from eqns. (3-10) and (3-12) 
using the same granulation trace to evalrmte the stored and real-time in- 
tensities, Ihe Standard MT (Ref. 6) was defined by: 


SN = 100, Al - IC^ 

^ rms ' 

M ^ = 66 (3-15) 

e± 

d = 0.05”. 

A large number of trackers wlxh variations on these standard values were 
studied three times each, using ground data, 65 cm, and 125 cm enhanced 
data. 

Figure 3-5 gives an approximate summary of these results, "while detailed 

plots are shoTm in Figures 3-6 xhrough 3-11. On all of these plots, errors 

are given in hundreths of an arcsecond. The solid line represents ground 

data (upper limits) and the dashed and dotted lines represent 65 cm and 

125 cm enhanced data, respectively (lower limits) . Noise errors for the 

standard MT lie between O.OI8" and 0.034", depending on the type of data, 

and systematic errors are between O.OO7" and O.OI9". Therefore, total 

tracking errors for fixed scenes are in the range O.OI9" to 0.039". 

Noise errors for the W} are relatively insensitive to SN, varying as 
-0 5 

SN . This happens becaus'e the quantization interval Al is much larger 
than the detector noise. Quantization acts as a nonlinear pre- amplifier. 
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Fig. ' 5-5 
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e GROUND DATA: UPPER LIMITS 



2 -SIGMA ERROR « 2 , UNLESS OTHERWISE STATED 



<9 ENHANCED DATA: LOWER LIMITS 



UPPER NUMBER FOR 65-CM DATA, LOWER FOR 125 CM 



SIGNAL-TO-NOISE RATIO 


Fig. 3-6 MT Noise Errors ( Centi-Arcseconds ) . From top to bottom, curves show results from 
ground data, 65 cm, and 125 cm enJianced data 
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Fig. 3-7 MT Noise Errors Versus Granular Contrast 
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Fig. 3-8 MT Noise Errors Versus N^j 
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Fig. 3-10 MT Systematic Errors Versus 


SYSTEMATIC ERROR 



converting small intensity changes to large^ fixed increments in Other 

simulations show that finer quantization^ say 10-bit instead of 3-bit^ would 
lead to larger tracker errors when SIT ~ 100 but improved performance if SET 
were Increased to 5P0. Figure 3-7 is included to show that a change in 
granular contrast_, AI , is equivalent to a change in SIT, only their product 
being relevant. 

-0 5 

Both noise and systematic errors decrease roughly as IT^^^ as shown 
in Figures 3-8 and 3-10; the slower decrease for ground data is a result of 
correlation peak broadening as IT^^ increases. The noisy behaviour of the 
systematic errors, which are averages over 100 independent scans, shows 
the large variation in tracker performance from scene to scene. This 
sensitive scene -dependence would also have appeared in Figures 3-9 and 3-11, 
showing dependence on image scale: however, these are smooth curves fitted 

throu^ three (five for ground data) simulation points. Recall that d, the 
detector elament size in arcseconds, is also a scale factor for the effective 
focal length. Errors are surprisingly insensitive to d over the range 
plotted, but increasing d must eventually cause very large errors '^dien 
elements become too large to resolve dark lanes and bright granules. 

(c) Time Dependence of Errors 

Two additional items of information are needed before the model of 
Section 3. 3 can be applied to the MT. First, the Standard MT scan pattern 
is a circle 1" in diameter. Therefore, althou^ the tracker length is 
3.14" or ~ 1,5 granule diameter, the MT only sees one independently moving 
feature: m = 1 in eqns. (3-7) and (3-8). Second, Appendix A shows that 

cr^ (At) = (o) e^^/'^o, (3-l6) 

where At is the time since the last update, cr^(o) is the fixed-scene noise 
error, and is the correlation time (Ref. l), 6,27 minutes. 

In all plots of time dependence, errors are in cent! -arc seconds and 
times are’ in minutes. Solid curves show upper and lower limits to absolute 
errors from eqns.’ (3-7) and (3-8); dashed curves show feature -following errors 
(eqn. 3-4), 

Errors for the Standard MT with groiond data are shown in Figure 3-12 
for times between zero and ten minutes. The update time, 30 seconds, is 
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Fig, 3-14 MT Errors for 125 cm Data 
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Fig, 3-15 Errors 
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Fig. 3-16 MT Errors After 1 Minute 




chosen to minimize the absolute tracking errors^ the solid lines. More 
informative are plots of the expected tracking errors after a fixed time 
interval (say^ ten minutes) as a function of update time as in Figure 
3-13. It is evident that feature -follo^ri.ng requires w 3 mins, for 
optimum performance^ whereas absolute tracking calls for in the range 
of 30 secs, to one minute. Regardless of the Standard M? shows large 
absolute tracking errors between 0.3*' and 0.7" after ten minutes. 

Improved telescope resolution^ using 125 cm enhanced data^ reduces 
feature-following errors but has no effect on the drift rates, as seen from 
Figure 3-14. The only way to obtain absolute tracking near the noise limit 
is to increase the tracker FOV to include many granules. Figure 3-15 plots 
the same error curves when the MT uses 66 elements 0.5" long around a circle 
of 10'* diameter (i.e., the Standard MT with d increased to 0.5**); I6 
independent granules are assumed in this FOV. Although the noise errors are 
increased, absolute tracking to better than I/4" over a ten minute period 
is possible. 

In some applications, such as use in conjunction with a second guider, 
the MT may only be required to track for short periods of time. Errors of 
the standard MT after one minute are plotted in Figure 3-16: updating is 

undesirable, and drift is only slightly larger than the fixed-scene errors. 
Even so, the tracking error is still comparable to a 125 cm telescope' s' 
resolution and is therefore unacceptable. 

3.5* Modified Ifershall Tracker (kMT) Simulation 
(a) lyfe-uhematical Description 

A simple change in xhe MT algorithm has been found which reduces the 
fixed-scene errors considerably. Specifically, the MMT uses the "squared 
residuals function" R(y' ) instead of xhe correlation function C(y’ ) used 
hy the MT. This function is defined by 

P 

E„ (y) = r [I. (y ± a) - (=y) ] , (3 -it) 

d 

where the symbols all have the same meanings as for the MT (see eqn. 3-8). 

It should be possible to convert the MT hardware to the WH algorithm with 
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minor changes in the electronics. Figure 3-3 shows examples of the residual 
function as measured hy the Standard MT (or MMT - the detector arrays are 
identical) and a CCD tracker described in Section 3 , 6 . 

The MMT locates a minimum in the residual function by the same method 
which the Ml uses to find -a correlation peak. Although noiseless, fixed 
scenes can be matched perfectly in principle by looking at the residual 
function, the MMT ^rill have a non-zero systematic error when R is not a 
symmetrical function of displacement. The MMT stable point is defined by 
(Mq) = systematic error is again given by eqn. (3-10). 

For fixed-scene noise errors, nhe MT equaxions (3-12) and (3-13) become 


2 2 (R) 

“ " 1 an/sy’ ' 


(3-18) 


o^(R) = 2 p AI^. (3-19) 

The quantities Al and p are again the quanxizaxion unit and the probability 

of error (see eqn. 3-lif). Equation (3-19) is derived in Appendix A and is 

true only for fixed scenes with residual functions much broader than the 

element size d. Jfoise errors become much larger as the two scenes and 

I.' decorrelate, as discussed below. 

*1 ^ 


(b) Fixed-Scene Results 

For computer simulation, the Standard MMP was assumed to use the same 
detector array as the Standard MT: SN = 5^0, = 66, d = 0.05”. 

summary of the fixed-scene errors is given in Figure 3-17* Noise errors 
for the Standard Mffl are smaller by a factor of 4 to 7 depending on the data 
than those of the MTj systematic errors show less dramatic improvement. 
Nevertheless, total errors for fixed-scene tracking lie between O.OO7" and 
0.010"; MMT errors are three to four times smaller for tracking over short 
time intervals (less than a minute) or for aligning two identical scenes for 
two different instruments. 

Figures 3-l8 through 3-22 give more accurate results than the rough 
formulae of Figure 3-17. Again, all plotted errors are in cenxi-arcseconds , 
and xhe three curves represent daxa-seta of different resolution. In general, 
the trends are similar to xhose discussed above for the MT. One curious 
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Pig. 3-18 MMT Noise Errors (Centi-Arcseconds). Results Prom Ground Data, 65 cm, and 125 cm 
enhanced data 
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Fig. 3-19 MMT Noise Errors Versus N^j^ 
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Fig. 3-21 MMT Systematic Errors Versus 



Image Scale 



difference is that Mffl! systematic errors are essentially the same for all 
types of data used; this accounts for the reduced uncertainty for the 
total fixed scene errors. 

(c) Time Dependence of Errors 

Feature-follo'wing and absolute tracking errors are again calculated 
from eqns. (3-4), (S'^Y) (3-8). For the Standard MMT, one feature in 

the FOV is assumed (m = l). As the real-time scene changes, perfect 
cancellation of the stored scene in the residual function no longer occurs, 
and noise-signal cross-terms cause xhe noise error (At) to increase 
rapidly. Appendix A shears that 

(At) « (o) e ^'*^'^'^0 [1 + 20.5 (1 - ]. (3-20) 

Figure 3-23 shows the f eatiire-foUowing and absolute errors (upper and 
lower limits) after ten minutes for the Standard MMT. Because of the rapid 
growth of noise errors, best update times are somewhat shorter than those 
for the MT. More frequent updating causes the feature-following errors 
and the most optimistic absolute errors to be about half those of the MT, 
although fixed-scene errors are three times smaller. The full advantage of 
the MMTis obtained for short time intervals, as comparison of Figures 3-24 
and '3-l6 shows. 

Drift rates are of course identical for xhe MMT and the MT. They are 
greatly reduced by increasing xhe FOV ten-fold, as shown in Figure 3-25. 

¥ith this image scale, absoluxe tracking to at least 0.2" accuracy over 
ten minutes is possible. For longer time, errors increase roughly propor- 
tional to the square root of the elapsed time. 

3.6. CCD Tracker Simulation 
(a) Jfetheraatical Description 

The CCD tracker which has been studied uses a detector array of two 
perpendicular lines of elements, which are analyzed independently. As 
usual, only one line was simulated and errors were increased by a factor 
of y? . Basically, xhe CCD tracker performs a least-square fit of the real- 
time scene to the stored scene, finding the position of minimum residuals. 
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Fig. 3-23 MMT Errors After 10 Minutes ( Feature -Poll owing Errors and Bounds 
for Absolute Error) 
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Fig. 3-25 Errors for MMT With FOV Increased Tenfold 



The residual functior. is defined, as in eqn. (3-17), 


R(y) = [!.’ iy, t> - I (0, 0) f . (3-21) 

i=l ^ . 

Tracker stable point is the point where ^(yQ) is an absolute minimum. 

For simulation purposes, the detailed method of varying y to search 

for uhe stable point has not been specified: error results are therefore 

somewhat idealised lower limits to an actual CCD tracker’s performance. 

However, it should be possible to approach these limits rather closely with 

an intelligently designed search algorithm and servo system. The real-time 

intensities I.’ (y, t) are measured at evenly spaced points y. = y + (i - l) d, 
i ^ 

and so the approximate location of the stable point can be found from each 

measurement of the , Note that this requires the computing power to 

evaluate the residual function several mimes, using data interpolated from 

the measured I.’ , for each real-uime look at the scene. 

1 

The CCD tracker has two conceptual advantages over the MT, in addition 
to larger values of SN and First, quantization effects can be ignored, 

assuming ten bits are available, which allows the tracker to take full 
advantage of its increased sensitivity. Second, it is clear that the sys- 
tematic error cr is always zero (recall that a is the tracking error on 
s s 

fixed, noiseless scenes). Only noise errors and tracker drift need be 
considered. 

The noise error for fixed-scene tracking is defined as follows. 

In the vicinity of y = y , the measured residual function R(y) is noisy with 
2 ^ 

variance a (R), v/hich is calculated is Appendix A, This noise causes nhe 
absolute minimum of R(y) to deviate from y . Then nhe rms noise error u 
is one-half the interval Ay such that 

R (y^ i Ay) - R (y^) = a (R), (3-22) 

where a (R) is given by eqn. (AI 7 ). Errors reported in the next section 
are ensemble averages of ct^, calculated from eqn. (3-22) using noise-free 
residual functions for 100 granular scenes. Use of one-half Ay as the rms 
error was justified by Mi3nte Carlo calculations: repeated addition of 

computer -generated noise to a variety of residual functions showed this to 
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Ise the correct one-sigma variation of the apparent stable point. 

Noise errors for tracking time-dependent scenes are defined in the same 
■way. To evaluate them, the Temporal behaviour of R (y, t) musT be knom. 

By expanding the squares in eqn. (3-21) and substi'tuting eqn. (3-3), the 
following approximate expression is found: 

R (y. t) [21 (o, o) + e“^/'"o I ' (y, o)] I, (o, o) (3-23) 

i IX 

In the vicinity of y = y^, the residual function has a parabolic shape for 
all trackers and data-sets studied: finite telescope resolution and detector 

element size always blur the intrinsic cusp-like behaviour expected of 
random granular scenes. Therefore, this result can be compactly rewritten as 

R (y, t) ^ R (y^, t) + l/2 e'^/^o r" (y _ y^ (3_24) 

■When this is used to e'valuate cr^ (t) according to eqn. (3-22), it turns 
out that 

(t) a e"/^o (R (t) ). '(3-25) 

. Finally, eqn. (AI7) of Appendix A shows that 

a/ tt) - a/ (o) [1 + 2 (SN - e"’^/'^°)] (3.-26) 

Noise errors increase rapidly in 'cime after each update, similar to the MMT 
behaviour. 

Tracker drift rate is calculated as usual from eqns. (3-7) and (3-8), 
providing upper and lower limits, respectively. 

(b ) Fixed-Scene Results' 

Fixed-scene noise errors were evaluated from eqn. (3-22) for a large 
grid of CCD trackers using ground, 65 cm, and 125 cm enhanced data. The 
Standard CCD tracker for illus-cration is 

SN = 500, = 1(4 

N^^ = 512 (3-27) 

d = 0.167". 
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Note that is the number of elements in each of the two lines of the 
detector array. Ihe detector element size d is the angiil.a.T size of a 
0.001 inch element in the focal plane of either a 125 cm F25 aperture or a 
65 cm F50 aperture; values appropriate to 125 cm F50 Q'lid 65 cm F25 were used 
in other calculations. 

Figure 3-26 gives a brief summary of the fixed-scene noise errors. 

OZhe Standard CCD tracker has a noise error between 0.0021" and 0.0004". 

The lower limits from enhanced data seem suspiciously small, even though 
they were computed in precisely the same \iay as the upper limits from ground 
data. In any case, these errors are between 16 and 40 'times smaller than 
■chose of the MT, and between 4 and 6 times less than the MMTls errors. As 
shown above by eqn. (3-26), they increase rapidly as the real-time scene 
deforms . 

The detailed plots shp-w noise errors in centi-arcsecond's versus SN, 

, and d for -the three data-sexs. As discussed above, the parabolic 
behaviour of R(y) causes very precise SN ^ dependence of the noise error. 
Quantization does not change this unless the interval Al becomes comparable 
to I /SN, which does not happen unless the number of levels is less xhan 

0<V 

about 5 SN*AI . In Figure 3-28, N ^ is varied between 20 and 220, a limit 
rms el 

set by ccanputer memory: the excrapolation to 512 elements seems safe. 

Variations of "with d in Figure 3-29 are also expecxed and unexciting. 

(c) Time Dependence of Errors 

The theory of feature-follo'wlng and absolute tracking errors is summarized 
in eqns. (3-4), (3-7) and (3-8). Featxire-folloTrlng errors are assumed to be 
caused by noise alone. Growth of the noise error o"j^(t) as the scene changes 
is given by e'qn- (3-26). Total length of the Standard CCD tracker is 85.5" 
in each arm, which was assumed to conxain 43 granules. 

Figures 3-30 and 3-31 show the absolute errors (upper and lower limixs 
in solid curves) and the feature follo'wing errors after 10 and 60 minutes 
obtained with ground data. Optimum update times are again between 30 and 
60 seconds for absolute accuracy and 2 to 4 minutes for minimizing noise errors 
alone. The worst-case errors are approximately 1 and 3 resolution elements 
after 10 and 60 minutes for a 125 cm aperture. Although many more curves 
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Fig. 3-28 CCD Noise Errox’s Versus 
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Pig. 3-29 CCD Noise Errors Versus Image Scale 
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Fig, 3-30 CCD Errors After 10 Minutes (Feature-Following Error and Bounds for Absolute Error) 
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could be generated, as for the MT and MMT, these two plots are sufficient 
to show the capabiliojr of the Standard CCD tracker. 


3.7. Forked Example of Tracker Error Analysis 

In this section, xhe errors of a hypiothetical CCD tracker are calculated 
in detail, to illustrate the use of the formulas and plots of the preceding 
sections. Although the example is completely contrived, none of the numbers 
chosen is unrealistic. 

Ff'oblem : A correlation tracker has two crossed lines of 256 detector 

elements each. The CCD array has elements 0.001 inch square to be used in 
the focal plane of a 125 cm F50 nelescope. The tracker performs a least 
square fit between its stored and real-time scenes, pointing the telescope 
to minimize the residuals of the fit. The real-time scene is measured 200 
mimes per second. If the detector elements have a signal-to-noise ratio of 
800, how accurately can it track granulation scenes at disk center? Speci- 
fically, calculate the rms tracking errors expected one and ten minutes after 
initial acquisition of the target. 

Solution : The first step is to calculate the rms tracking error caused 

by noise on scenes which are independent of time. Upper and lower limits 
for this quantity, are most easily obtained from the equation in Figure 
3-26. Substituting SN = 8OO, Al^^jg = granulation contrast ax disk center = 
IC^, = 2 %, we obtain 

= 0.77 ( 037"^°'^^ '"o < .0021". 

Here, d is the angle subtended by the length of one detector element in the 
focal plan: 


d 


, 0.001 inch \ 
^ focal length' 


O.OSi!- arc seconds. 


This leads to 0.0002" < < 0.0012"; the upper limit will be used in the 

rest of the calculation. Errors in both x and y directions are included 
in 

If the tracker were always used to align two identical scenes, then 
the only trackX 3 :ig error woiild be ^ 0,0012". However, this result could 


5-^7 



he changed by the frequency response of the servo loop of which the tracker 
is a part as discussed in Section 2.7* As an error sensor, the tracker has 
a bandvridth of 0-200 Hz with a t/hite noise spectrum. However, since the 
tracker uses crossed lines instead of a raster pattern, it would lose its 
lock on the target if image displacements of more than 1" or speeds 
exceeding 200 arcsec/sec were encountered. 

In most applications, the tracker would be attempting to match time- 
variable scenes. Changes in the noise error in time and the tracker drift 
rate caused by motion of granules in its ¥0Y must be considered. The noise 
error cr^(t) is given by eqn. (3-26), where = 6.27 minutes is the granular 
correlation time, Ibr t = 1 minute, this equation gives ct (l) = 0.0064"; 

n 

this is the "feature-follov?ing" error discussed in previous sections. 

The absolute error, includiixg noise and tracker drift, is bounded by the 
results of eqns. (3-7) and (3-8). Substituting V = rms granular speed 
= 1 km/sec - 0.083 arcsec/min, m= number of granules in FOV - N ^*d/2" = 

10, and n = number of updates = 0 , we find that the tracker will drift by 
about 0,026" in one minute, on the average. This is much greater than the 
noise error, and so the total absolute error afuer one minute is 

E (l) 0.026" rms. 

To track accurately over uen minutes, the tracker must update its snored 

picture periodically. For absolute tracking, an update time t =1 minute 

is near the optimum value, as shown by Figure 3-30. The accumulation of noise 

errors is now described by eqn. (3-4), with a = o, n = 10, and At = o; 

s 

this yields 

(10) = 10 (0.0064)^ + (0.0012)^; 

E^ (10) = 0.020”. 

Upper and lower limits to the total error come from substitution into eqns. 

(3-7) and (3-8): 

0.085" < E (10) < 0.21". 

Thus, the tracker will drift by one or two resolution elements in a ten 
minute interval. 
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Two final comments are in order. Firsts = 1 minute is the average 
update time. In fact^ the tracker should update when its Tnimmiim residual . 
rises to a fraction 

(1 - e"'^u/To) 0.15 

of the uncorrelated value;, 2 <lS’. Second the absolute tracking errors 

refer to coordinates moving with the average solar rotation at disk center_, 
0.15 areseconds/inin. If desired, the rotation could he removed independently 
of correlation tracker operation. 

3.8, Summary of Tracker Performance 

The three xrackers analysed in detail above are suimnarized in Figure 
3-32. This shows upper and lower limits to the rms tracking errors uen 
minutes after initial target acquisixion. The trackers are used tri.th a 
65-125 cm telescope to follow specific features or to track in absolute 
coordinates with respect to the solar limit. 

The present Marshall Tracker shows errors between one and ten telescope 
resolution elements. It suffers from a small field-of-view and a poor 
choice of algorithm, aggravated hy low sensitivity and three-bit quantization. 
The modified Marshall Tracker uses the same hardfra,re with a minor (but 
effective) change in algorithm and with a variable image scale. Its per- 
formance is quite good, with errors of at most two resolution eleiaents 
after ten minutes. Finally, the CCD Tracker corrects all shortcoming of the 
Itershall Tracker: sensitivity, FOV, algorithm and computing capability 

are all greatly improved. Taken together, these improvements result in 
errors one-third to one-half of those for the Modified Marshall Tracker. 
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TRACKER PERFORMANCE SUMMARY 




BLOCKBIEe.O 




ERRORS (ARCSEC) AFTER 10 MIN 

TRACKER 

FEATURE-FOLLOWING COORDS. 
UPDATE TIME 1-2 MIN 

ABSOLUTE COORDS. 
UPDATE TIME « 30 SEC 

MARSHALL 

0.07 < € <0.10 

0.20 < e < 0.70 

MODIFIED 
MARSHALL WITH 
VARIABLE IMAGE 
SIZE 

0.02 < € < 0.07 

0.06 < e< 0.20 

STANDARD CCD 

0.007 < e < 0.03 

0.02 < e < 0.10 
















4. Focus Quality Informtion 

4.1. Introduction 

In principle^ the detector array of a correlation tracker can he used to 
measure the focus quality of its image. Defocus, misalignment or aberrations 
al'ways decrease the fine structure or high spatial frequency power in the 
image. Suitable processing of hie tracker detector outputs can measure this 
fine structure content. Biere are tiro potential uses for such a system: 

(1) to assist in finding the position of sharpest focus, during a focussing 
operation; (2) to act as a constant monitor of image quality, signaling 
when refocussing is necessary. For either application, a single parameter 
derived from the tracker scene is needed for use in a control loop. The 
scene-to-scene variation of any parameter, even when focus is perfect, 
presents a problem for the second application; the optics must be focussed 
on each scene, and refocussed periodically as the scene changes, to provide 
a baseline for sensing image degradation. 

In this study, two different indicators of focus quality were evaluated 
using statistical descriptions of granulation scenes. Detector configurations 
corresponding to the MI and the CCD tracker were modelled, and the average 
sensitivity to noise and to defocus of hoth indicators were computed. These 
indicators are: (l) the curvature of xhe correlation function, which can he 
obtained from the present MT id.th trivial modifications; (2) the image 
power in a pre-selected spatial frequency hand, which requires real-time 
calculation of Fourier transformed images. This s'econd indicator is an 
effective measure of image quality if used with a CCD tracker, but it is 
not a panacea for the problem of active focus control. 

4.2. Jfe,thematicaZL Ifodel 

The mathematical formulation used to study focus quality indicators is 
described in this section. The fundamental theory is first presented as if 
a continuously sampled image were available, for simplicity. Then the effects 
of the finite number of detector elements and their spacing are considered. 

All calculations are one-dimensional, as if detector arrays were straight 
lines: this is a good approximation as long as the spatial frequencies 

considered are longer than the inverse curvature of the tracker element; pattern. 




Given a one-dimensional image l(x) ^ i;he MT measinres an approximation 
to the correlation function C(y) defined by 

1 V2 

"" L J (^-1) 

-L/2 

’where L is the total length (in arcseconds) of the detector array. !Ehe cujrv- 
ature D tjhich can be measured by the MT is 

D = 1 C(d) + C(-d) - 2 C(o) \., (k-2) 

d being the element size. Defocus always causes the cu 2 Tvature to decrease. 

If a tracker uses the residual function B(y)^ it is easy to show that the 
curvature is the same, on the average. 

The second image quality indicator is based on the Fourier transform 
of I(x): ^ 

I (k) = J g2rrikx ^^^2) 

-CO 

here, k is nhe spatial frequency with units cycles/arcsecond. Power in a 
band, k^ < k < k^, is defined by 

P (k^, ^2^ " J 1 I (fe) dk. (Vi)-) 

By proper choice of k^ and k^, power in a band can be a sensitive measure 
of image quality. The value of the correlation peak curvature can also be 
expressed in terms of I (k) by means of the relation: 

03 ^ 

D o: J [1 - cos 2 TT kd] | I (k) dk^ (^“5) 

—CO 

the exact proportionality constant will not be needed. 

Actual (as opposed to measured) intensities of a granulation scene, 

(x) or (k), can be described statistically using the power spectrum 
(k) . Specifically, 

< I (k) |^> a (k). (4-6) 
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Sq (k). as measured by Kinahan (Reference 2, Chapter 3}^ is plotted in 
Figure 3-4. She intensities measured by a tracker detector array are described 
by the modified power spectrum S (k), 'which includes the effects of telescope 
resolution and detector element size: 

S (k) = (k) j MiF (k) T (k) |^. (4-?) 

MTF (k) is the modulation transfer function of the optical system, including 

r>^ *■ 

defocus or aberrations, and T (k) is the Fourier transfom of the detector 
element profile. 

i („) = <' . (i»-8) 

' tt k d 

All of these functions are even funcoions of k, and so only positive spatial 
frequencies need be considered, if desired. 

It is convenient to express MEF (k) as a product of two terms, 'where 
the first represents the perfect, diffraction limited MEF, and the second 
simulates loss of resolution caused by defocus. Ihis is nox a rigorous 
procedure, by any means, but an exact treatment would require a much more 
elaborate specification of the optical system and its aberrations. Ihe 
follo-tang simple model is certainly adequate for a one-dimensional calculation: 

MTF (k, ct) = MTF^ (k) • A (k, o) , (^-9) 

where 

(k) = 1 - (4-10) 

A (k, cr) = exp [- TT (ct k fy . (4-U) 

The largest detectable spatial frequency is 

k = (4.85 cy/ ") (telecsope diamter/50 cm). (4-12) 

in3«^ 

Defocus is measured by the dimensionless pararaexer cr, which is defined by 
eqn. (4-ll): this definition was chosen so that the effective resolution 

is given by 
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, Ef A active \ Jt , ? 

( „ ^ . ) fa VI + cr 

' Resolution 


, Diffractions 
^ Limit ^ ' 


(4-13) 


Equations (4-7) through (4-l3) permit the calculation of < 1 1 (k) j^> 
for any values of the aperture size^ defocus parameter and detector 
element size d. The focus qualixy indicators can then be obtained from 
eqns. (4-4) and (4-5). Finite sampling effects are easily taken into account. 
If xhe tracker has elements spaced at intervals of Ax, then I (k) can 


only be measured at the discrete frequencies k^, xtiere 




k. = 

a 


lieiAx 


0 = - 


el 
2 " 


5T 


el 


(4-14) 


Each Fourier component I (k.) is an average over an interval Ak=l/l Ax. 

Finally, the sensitivity to detector noise of measured values of D and 

P (k^, kg) must be calculated. Let n. be the noise in the output of -the 
th ^ 

j — detector element. Since xhe various n. are independent, the noise power 

3 

spectrum has a consxant value, P^, over the allowed frequency band defined 


by the k.. 

3 


Using the fact that 
2 


< n“ . > 
3 


a 


S Ak Pjj (kj. 


(4-15) 


one can derive expressions for in terms of (k) , maintaining consistent 
units. The results are 


S S (k.) 

P = 

U 


U,(SNAI ) 
el rms 


( 4 - 16 ) 


for non-quanxized data, and 




0,44 S S (k.) 

,1 ^ ^ 


( 4 - 17 ) 


for three-bit quantized data. SN is the detector signal -to-noise ratio 

before quantization, and Al is defined in eqn. (3-2). One P„ is known, 

■* ,rms 2 a \ / jj f 

the standard deviation of | I (k . ) j is given by 

J 


g ( lir) 
li"l 




B(h.) 


S(k.) 1/2 

[2 + 4^ ] . 
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( 4 - 18 ) 



Noise sensitivities of D and P k^) are easily- computed from this formula. 


4.3. Results 

Correlation peak curvature D and power in a selected bandwidth P (k^^ k^) 
were calculated for two aperture sizes and two detector arrays. These 
results are most clearly summarized by a set of plots, shotring the image 
quality indicator vs. effective resolution. The behaviour of each indicator 

^ ' f 

in the vicinity of perfect focus is desired, and so curves' 'al'e' normalized 
to their dzffra'ct'ibh I'imited value's. Soil'd curves show expectation value's' , 
and dashed curves' are ± one-si^a variations to show the effects of detector 
noise. 

Figure 4-1 shows D (defined by eqn. 4-2) for the Standard J&Lrshall 
Tracker of Section 3.4 used with a 50 cm aperature. The decreasing curvature 
with increased blurring is evident, but noise renders this useless as a focus 
indicator. A Standard CCD array (Section 3.6), as plotted in Figure 4-2, 
measures D accurately, but a sharper peak at perfect focus would be desirable. 
Although simple to compute and atxractive in principle, D does not appear, to 
be useful in practice-. 

Pbiirer in the 1-2 cycle/arcs econd band is shown in Figures 4-3 and 4-4 
for the two trackers. This bandwidth was chosen xo show xhe best possible 
behaviour of the Marshall Tracker as a focus quality sensor. Extending the 
band xo higher frequencies admits too much noise, whereas a lower frequency 
band is less sensitive to focus changes. The figure shows that noise errors 
are comparable to a 5C^ loss of resolution: without a sizable increase in 

SN, xhe Ifershall Tracker is a poor instrument for measuring image quality. 

On the other hand, a CCD tracker has the sensitivixy and xhe spatial 
resolution to make use of even higher spatial frequency components, if desired. 
Figures 4-4 through 4-6 demonstrate the potential for very accurate location 
of sharpest focus by adding Fourier transform capability to a CCD correlation 
tracker. The curve in the last figure is less steep than the others because 
detec-fcor sizes are nearly unchanged while diffraction limited resolution has 
increased by a factor of 2.5* Even so, the tracker can still sense resolution 
losses of 5^ or less. 
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Conclusions and Recommendations 

This study has characterized correlation trackers in three distinct 
i/ays. Chapter 2 gives an extensive^ qualitative discussion of the utility^ 
design principles^ key parameters^ and error sources of correlation trackers 
as applied to pseudo-random images in general and the solar granulation in 
particular. Three specific tracker concepts are analyzed quatitatively in 
Chapter using observational data on the solar granulation to calculate 
expected tracking errors in various potential applications. Pinally^ 
chapter 4 sketches the ability of a correlation tracker to perform a 
second independent task, the direct measurement of image quality at the 
focal plane. 

In general, a correlation tracker is appropriate for image motion 
compensation when the images are random, only their statistical properties 
being knoim in advance. The solar granulation is a random, low contrast 
image with typical correlation lengths and times of one arc second (725 tm 
on the sun) and six minutes, respectively. Its temporal statistics are 
not well known, and so analyses of the time-dependence of tracker errors 
are necessarily tentative. Nevertheless, a well-designed correlauion 
tracker can perform acceptable image motion compensation on a meter-class 
solar telescope in three possible configurations : 

a Use within the loop of a limb guider; the tracker stores images 
whenever the guider error signal is small and performs active 
tracking only for the short -Dime inxervals when uhe guider error 
signal exceeds a threshold. 

e Alignment of identical images fed to two different instrumenxs; 

« Tracking in absolute coordinates wiuh respect to the llmbj 

although a CCD tracker can have errors comparable to telescope 
resolution for tens of minutes, it should not be used as the only 
pointing system. It cannot acquire a target Initially, reacquire 
after a displacements greater than one arcsecond, or remove the 
effect of solar rotation. 

Three sources of correlation tracking error have been identified and 
studied. Errors caused by detector noise can be minimized by proper choice 
of the tracker algorithm, number of quantization levels, detector sensitiv- 
ity, and number of detector elements. Systemtic errors are image dis- 
placements which occur when aligning tiro identical, noise-free scenes; 
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they can he eliminated by choosing a proper algorithm or minimized by 
increasing the number of elements. Tracker drift, caused by the motion 
of granules in the FOV, can only be minimized by increasing the FOV and 
choosing a proper scan pattern. 

The numerical simulations of chapter 5 are summarized in Figures 
5-5^ 3 - 17 ^ "and 5-26 for tracking fixed, time -independent scenes; the 
growth in time of tracker errors is summarized in Figure 3-32. The 
present Marshall Tracker (Ml) does fairly well on fixed-scenes but drifts 
very rapidly because it only sees one granule in its FOV. If the M? is 
modified to use the least-squares residual algorithm, its fixed scene per- 
formance improves significantly. Additional inclusion of a variable image 
scale can easily cut the drift rate by a factor of three. Finally, the 
standard CCD tracker shows very small errors for all applications listed 
above , 

Chapter 4 sho^/’s that a correlation tracker can be used to measure 
the focus quality of the image, if Fourier transform capability is added 
and if the nracker scan pattern includes a uniformly sampled linear array. 
The Jfershall Tracker is of doubtful utility for this purpose because of 
its low detector sensitivity and three-bit quan-cization. A CCD tracker 
can measure focus quality very accurately, as shown in Figmres k-k through 
4-6. However, the following reservations must be noted: 

9 the tracker cannot tell which degree-of -freedom of the optical 
system is causing the defocus; 

e the tracker cannot always distinguish between image scale changes 
and defocus; 

0 the tracker's image quality indicator is scene-dependent, and so 
frequent updating may be needed during refocussing; the interaction 
of refocussing and simultaneous tracking deserves more study. 

There are several areas of investigation which merit further study. 
First, since CCD trackers exhibit the greatest promise in all performance 
areas except scan pattern versatility and low light level tracking, a 
detailed design of a CCD tracker should be undertaken, and a tracker con- 
structed. Such a tracker should also be designed to serve as a focus 
monitor. 
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Second^ since the detailed characteristics of the time varying 
granulation patterns are at present unknovn^ this study has left rather 
wide error "bars in the errors ei^pected as a function of time. Results 
recently obtained by Dr. P. Mehltretter could be used to vastly improve 
these results^ and could be incorporated in a follow-on study. 

It is perhaps wort^h pointing out that an engineering design of a 
correlation tracker specifically for use with the SOT would be useful for 
producing heat source, weighu, and size numbers for SOT engineering 
definition studies, even if a tracker is not actually constructed. Also 
the characteristics of the tracker as a servo component would be better 
defined. 

A study is also needed on interactions between the limb tracker and 
correlation tracker, although this will probably have to wait tmtil the 
telescope structure is better defined. A part of this problem which 
could ’be immediately attacked is the question of techniques for compensating 
for solar rotation. 

¥e recommend that development of image dissector trackers proceed. 

We do not believe, however, that the specific configuration of the 
Marshall tracker studied represents an immediately useful configuration 
for SOT. 
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Appendix A. Tracker Moise Errors 


1. Intro due tion_ 

This Appendix gives brief derivations of the foraiulas used in Chapter 
3 for the noise sensitivity of the three tracker algorithms. The important 
quantities are the variances due to noise of the correlation function C(y) 
and the residual function R(y), >ri.th either quantized or unquantized intensities. 
Detector noise level is assumed to be constant^ with standard deviation 

although shot noise may be dominant for some detectors, it is neg- 
ligibly different for a low contrast target like the solar granulation. 
Throughout this Appendix, all averages are taken over the possible values 
of the noise, and not over an ensemble of scenes 

2. Noisy Correlation Functions 

The detector outputs are two sets of intensities T^/ith noise, I^ + n^ 

and I.’ + n' . j these are quantized to eight levels after I subtraction. 

1 X 2 - ‘ av 

The variance o (C) is needed, where 

C = 2 (I + n^) (I’ + n^) (Al) 

and 

(C) = < C^> - < C >^, (A2) 

The noises n^ and n^ are independenx and have zero mean, and so 

<0=2 I. 1\. (A3) 

XX 

X 

After some algebra, one obtains 

< C^ > = 2 I. I . I! T. (A4) 

i,- X a X 0 

+ < 2 nj r.^ + + n? nl%>. 

. i X XX XX 
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When the intensities are quantized into intervals Al much larger than the 
detector noise a good approximation is 

< n ^ > = p AI^ = ( — ^ ) AI^. (A5) 

SN AI 

p is the probahility that the measurement ■tra.s mis-quantized because of 
the noise n^, and p is a numerical constant of order unity. Assuming the 
intensities are uniformly distributed across each quantization interval, 
a simple calculation using the properties of the normal error function 
yields the value O.698 for p. 

Assembling all of xhese equations, the variance becomes 

(C) = p Al^ 2 (I,^ + r.^ + P AI^), (A6) 

■which is essentially the same as eqn. (3-13). The degree of correlation 
between and is clearly irrelevant, and so this result is not limited 
to fixed scenes. When the scenes decorrelate in time, the change in 
is caused by the exponential decrease in correlation peak slope, and so 

a (t) = a (0 ) (At) 

3. Woisy Residual Functions 

Using the same notation as above, the residual function R is given by 
E = 2 (I. + n - l! - n’)^. (A8) 

i 1 1 3 . 1 

Again, n. and nl are assumed to be independent wath zero mean. I. and l! 
may be identical (in the CCD algoriubm with fixed scenes), or they may be 
only partly correlated. In general, let p be defined by the relation 

S (I. - l’.)^ = 2p S l/. (A9) 

Then the general expression for the variance of R can be found by tedious 
algebra: 
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(AlO) 


(R) = S 8P < (n^ - n^. > 


/ _t\2 2 

^ • 


+ < (n^ - n^) > - < (n^ - 


In the MMT algorithm, all intensities are quantized; therefore, the 
noise averages are evaluated using eqn. (A 5 ) and 


< n^^ > = p 


(All) 


These lead to the result 


O^ (R) = 2 p Al^ + 2 p Al^ 


+ 16 p p AI^ S 
i 


'(A12) 


When the MMT operates on fixed scenes, p 0.02, and the first term of eqn. 
(A12) is dominant. The second term is always negligible, hut the third 
grows rapidly because, from eqn. (3-3), 


p 1 - 



(A13) 


Using the rule-of-thumb that the rms fluctuation is approximately one-fifth 
of the peak-to-peak, it folk)WS that 

S s ( I AI^. (A14) 

i 

Assembling eqns. (A12) - (A14), the final result is 

(R) = 2 P AI^ C It 20.5 (1 - ]. (AI 5 ) 

This accounts for eqns. (3-19) and (3-20). 

Quantization effects are ignored for the CCD tracker, and normal 
distributions are assumed for the noise. In this case, 

< n^^ > = 3 < >^ = 3 ( ^ )^. (AI 6 ) 

Combining this with eqns. (AlO) and (A13), the variance can be written as 
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/ (E) . 8 ( ^ 


el " SW 


2(1- 2 I,^ ], 


which justifies the treatment of Section 3,6 (a). 
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